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Abstract 

Hadron spectra and other properties of quark systems are studied in the framework of a 
non-relativistic spin-independent phenomenological model. The chosen confining potential is 
harmonic, which allowed us to obtain analytical solutions for both meson and baryon (of equal 
constituent quarks) spectra. The introduced parameters are fixed from the low-lying levels of 
heavy quark mesons. The requirement of flavor independence is imposed, and it restricts the 
possible choices of inter-quark potentials. The hyper-spherical coordinates are considered for the 
solution of the three-body problem. 

Espectros de Hadrons e outras propriedades de sistemas de quarks sao estudados do ponto de 
vista de um modelo fenomenologico nao-relativi'stico independente de spin. O potencial confinante 
escolhido e harmonico, o qual nos permite obter solugoes analfticas tanto para os espectros dos 
mesons quanto para os dos barions (de iguais quarks constituintes) . Os parametros introduzidos 
sao fixos a partir dos primeiros estados excitados dos mesons pesados. A condigao de independencia 
de sabor e imposta, o que restringe as possfveis escolhas de potenciais inter-quarks. As coordenadas 
hiper-esfericas sao consideradas para a solugao do problema de tres corpos. 

PACS numbers: 03.65G, 14.40, 14.20 
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I. INTRODUCTION 

The study of the fundamental or constituent blocks of matter has been for long time a 
fascinating field in physics. With the pass of the years new fundamental blocks or particles 
have appeared, modifying old concepts. For example the atom, that was initially supposed 
to be fundamental, was found to be formed by the nucleus and electrons, the nucleus was 
found to be formed by neutrons and protons (nucleons) and finally the nucleons to be formed 
by quarks. 

Nowadays it is believed that the theory describing the interactions between fundamental 
particles is the Standard Model. This model has a SU (3) x SU (2) x U (1) symmetry. The 
SU (3) component is what is called Quantum Chromodynamics (QCD) and it is the gauge 
field theory describing the strong interactions of quarks and gluons. The SU (2) x U (1) 
component is called Standard Electroweak Model describing interactions between quarks 
and leptons through the W bosons and the photons. There is an additional particle in the 
electroweak model, the neutral Higgs scalar that appears after the spontaneously symmetry 
breaking mechanism, whose existence has not been probed experimentally and that is a 
subject of a lot of present experimental research because it could prove the validity of the 
theory. 

We will be interested in the present paper in the low energy region of the QCD theory, 
in which quarks interact strongly to form bound states known as hadrons. When these 
bound states are formed by a quark and an antiquark (qq) they are called mesons; and 
when they are formed by 3-quark states (qqq) they are called baryons. The study of the 
hadron spectra is a fundamental and open field in theoretical physics. Up to the moment 
the more important approaches to this subject have been: Lattice QCD, QCD sum rules 
and potential models. Lattice QCD is the most fundamental approach and together with 
QCD sum rules have had good success. Potential models, although less fundamentals, have 
proved to be very useful even in non-relativistic approximations. Since the latest seventies 
a lot of attempts have been made in this field with very good results 

An important question for the use of potential models is whether it is possible to consider 
hadrons as non-relativistic bound states of quarks. This question could be answered solving 
the Schrodinger equation assuming a qq potential and then checking if the obtained quark 
velocities are non-relativistic. For the charmonium system this was done in Ref. and it 
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was obtained a result of (v 2 /c 2 ) = 0.2 in the ground state. For lighter hadrons the results 
are no so good and the validity of the non-relativistic approximation depends strongly from 
the interaction potential chosen. For example, the interaction potential proposed by De 

n it 

Ruiula, Georgi and Glashow [3j was shown 41 to be unsuitable for dynamical non-relativistic 

nnn 

calculations of light hadrons. However, Martin showed latter [5J, |6J, 17[ that a non-relativistic 
model with a power law potential is able to describe heavy quark mesons and the clearly 
relativistic ss states. The potential proposed by Martin had the form 

V (r) = A + Br a , (1) 

where A = -8.093 GeV, B = 6.898 GeV and a = 0.1. 

The study of baryon systems with this power law potential was done by Richard 
obtaining good results. At that time the baryon spectra in potential models, although 
rather elaborated jjj, was completely disconnected from the meson sector. Exceptions were 
some attempts to derive meson and baryon potential energies from a common framework as 
the instanton, string or bag model. The rule adopted by Richard for the qq potential was 

V m = \v m . (2) 

Other very well known phenomenological non-relativistic potential models have been 
based in the Cornel potential 

V (r) = A c r - ^ + C c , (3) 
r 

where A c = 0.1756 GeV 2 , B c = 0.52 and C c = -0.8578 GeV; and the logarithmic potential 
Ell 

V (r) = A L + B L In (r) , (4) 

where A L = -0.6631 GeV and B L = 0.733 GeV. 

Finally, other two propositions of inter-quark potentials are due to Song and Ling [1J 

V (r) = A s r^ + B s r~^, (5) 
where A s = 0.511 GeVi and B s = 0.923 GeV^- and Turin Q 

V (r) = -A T r"f + B T A + C T . (6) 
where A T = 0.620 GeV*, B T = 0.304 GeV* and C T = -0.823 GeV. 



Studies with phenomenological potentials, like the Cornell one, but considering a rela- 
tivistic kinetic ener gy term, are able to describe the observed spectra of heavy and light 
hadron systems 0,13] • Non-trivial connections between these relativistic potential models 
and rigorous numerical results from lattice QCD have been demonstrated |l6j . Some recent 
works have tried to understand why the non-relativistic treatment works and allows useful 

predl c tl o lls eve lltorrel a tl v lstlcsystems 0ia. 

In the previously mentioned potential models, central conditions have been the flavor 
independence of the potential chosen and the existence of a confining term. Quark masses 
appearing in these phenomenological models are the so called constituent quark masses 
that should not be confused with the current quark masses, that are the mass parameters 
appearing in the QCD Lagrangian. Constituent quark masses are bigger than current quark 
masses and it is suppose that this is due to gluonic condensate effects. In general constituent 
quark masses are considered as free parameters to fit in potential models, that is why different 
values are found all over the literature. Up to our knowledge there are two works in which 
constituent quark masses are calculated from QCD: one is due to Elias et al. |3| 
that used an operator-product expansion (OPE) of approximate non-perturbative vacuum 



expectation values in the fixed-point gauge; the other work is due to Cabo and Rigol 



2l| 



in which a Modified Perturbative QCD expansion incorporating gluon condensation was 



employed 



22|. Because there are no free quarks, a lot of care should be taken with the 



meaning of these quark masses. 

In the present work we study the hadron spectra within a non-relativistic spin- 
independent phenomenological model, with a harmonic confining potential. The idea is 
to show with an educational perspective what can be done with the Schrodinger equation, 
the very useful for physicists harmonic potential, and the same number of parameters used 
in almost all the calculations, in order to understand the hadron spectroscopy. The har- 
monic potential has the great advantage that allowed us to obtain analytical solutions for 
both the meson and baryon (with equal constituent quarks) spectra. Although the obtained 
results are not as good as the ones obtained for the former potentials, we think that they 
give very good estimates of the hadron properties without the need of numerical calcula- 
tions implemented for all the previously mentioned potentials. As we will work with a spin 
independent model, the magnitudes we will deal with will be spin averaged. For fitting 
our parameters spin averages of experimental values were calculated, but for certain reso- 
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nances the experimental values were not at hand [23( and theoretical results obtained by 
Fulcher 2^1 were used. It should be mentioned that in a work by Hirata et al. j^j] the 
harmonic oscillator was employed as an unperturbed confinement potential, in the asymp- 
totically free colored-quark-gluon model in which the one-gluon-exchange force was treated 
perturbatively. Later Ram and Hasala j3| used the pure harmonic oscillator potential in 
the Klein-Gordon equation to determine the meson masses. 

The exposition will be organized as follows. In Section |H] meson spectra is calculated 
for two different potentials with harmonic confining terms, that allowed to obtain analytical 
solutions. The first potential is the harmonic oscillator and the second one is the harmonic 
oscillator plus a term proportional to \. The results are analyzed for the flavor dependent 
and independent cases. In Section IIIII we study the properties of the radial wave function, 
and some related physical quantities, at the origin. Section IIVI is devoted to the study of 
the baryon spectra and the summary can be found in Section 5. Finally, an Appendix was 
introduced for mathematical details of the three-body problem solution. 

II. MESON SPECTRA 

In the present section, the non-relativistic meson spectra are calculated for two different 
potentials. The first one is the pure harmonic potential and in the second one a term 
proportional to \ is added to the harmonic oscillator, allowing to improve the short-range 
interaction. That is the potentials considered are: 

Jc-i V 2 

Vi (r) = ^- + W U (7) 



and 

k 2 r 2 a 
~2 r 



V 2 (r) = -L--^ + W 2 . (8) 



The Schrodinger equation, in the Center-of-Mass (CM) system and in spherical coordi- 
nates, has the form (notation h — c — 1 is considered) 



i r i a f r2 o_\ i 



2/x \r 2 dr \ dr J r 2 



1 8 ( n B\ 1 d 2 
sins'— + 



sin 9 89 \ 89 J sin 



,2 



{Ei-Vi)}% {r,9,(f>) = 0, 

(9) 



where i — 1,2 and \x = is the reduced mass. Introducing the spherical harmonics 

^i,n k rr H {r, 9, 0) = ^R i>nl . (r) Y i>km . (9, 0) , (10) 



the radial Schrodinger equation can be written as 



f 1 d 2 ky L (L + 1) ,) 

where 

1(1+1) for i = 1, 
I Z (Z + 1) - 2/ja for i = 2. 
The solution of Eq. fjX X j) can be found in any classical textbook of Quantum Mechanics 
(27, 28] and has the form 

*wo=jw=*W (-«,<, + §,«•), (12) 

where 

r 



1 



A^j jn ;. is a normalization factor and F (— n, Zj + |, is the confluent hyper-geometric func- 
tion. The eigenvalues of the energy are 




E i = J-\2n + l i + -J+W i . (13) 
First, in order to compare with experimental values we fitted the parameters without 



a 



it is only 



imposing the flavor independence condition on the potentials. For that case 
possible to determine from the experimental spectra the values of e — yf^, 2/xa and Vq = 
Wo + mi + rri2 that are shown in Table [H for charmonium (cc) and upsilon (bb) systems. As 
it can be seen, the parameter 2fia has the limiting value for the existence of the solution for 
the flavor dependent (FD) potential V2 for both quark systems, because 




and for I = and 2/xa > 0.25 the squared root will have imaginary values. 

The mass values obtained for cc and bb systems, with the parameters in Table HI are 
presented in Figure ^ together with the theoretical values obtained by Fulcher j3| and com- 
pared with the experimental values [3| (the results for each potential are presented in the 
figure as columns). The experimental masses of the states IS*, IP, 2S for the charmonium 
systems and all the ones presented for the Upsilom systems are spin averaged and were used 
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TABLE I: Parameters obtained for the cc and bb systems with the flavor dependent potential. 



System 


e (MeV) 


2fia 


V = W + m 1 + m 2 (MeV) 


Vf Da (cc) 
V 2 FDb (cc) 
Vf D (bb) 
V 2 FD (bb) 


296 
302 
217 
219 


0.25 
0.25 


2679 
2773 
9235 
9301 



"Flavor dependent potential V\ (r) 
h Flavor dependent potential V2 (r) 



4400 


3S- 


4307 


4283 


4200 






4106 


4000 


2P 


4011 




3800 


"1D 

2S 


3715 


3814 
3679 


3600 






3502 


3400 


1P 


3419 




3200 


1S- 


3123 


3075 


3000 




Vf D 


v 2 FD 



4252 



3966 



4081 



3915 



4040 



3763 3770 



3661 3663 



3397 



3467 



3525 



3068 



3112 3102 

Vf FU[29] Expt[23] 
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FIG. 1: Charmonium (left) and Upsilom (right) energies (MeV) (V FI is the flavor independent 
potential V\ (r)). 

in our fit of the parameters. The experimental masses reported for the ID and 3S charmo- 
nium resonances belong to the l 3 Di and 3 3 Si states. The results obtained for the potential 
V-i are in better agreement with the experimental results because the non-harmonic term 



t-2 , 



improves the behavior of the potential in the region r — > 0, and breaks some degen- 



eracy present in the solution of the pure harmonic oscillator. Unfortunately the potential V 2 
was found incompatible with the flavor independence condition, been impossible to obtain 
a unique potential like this for charmonium and upsilon systems. Then, after requiring the 
flavor independence condition only the pure harmonic oscillator parameters can be fitted, 
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and the following values were obtained 

k = 0.155 Gel/ 3 , m s = 2.725 GeV, 
Wo = -4.94 GeV, m c = 3.812 GeV, 
m b = 7.093 GeV, 

where the s quark was included as in Martin's works. The spectra calculated with these 
parameters, for cc and bb mesons are presented in Figure ^ (Vf 1 ). Other meson resonances 
like ss, cs — sc, bs and be are shown in Table ITT1 The mass of the IS" resonance for the ss 
meson was employed to fit the s quark mass. 



TABLE II: Other meson resonances (MeV). 



State 


Vf^ss) 


Expt a (ss) 




Expt a (cs) 


Vf^sb) 


Expt a (s6) 


Vf^bc) 


Expt a (6c) 


FV b {bc) 


IS 


1016 


1019 


2065 


2076 


5299 


5370 


6340 


6400 


6361 


IP 


1353 




2377 




5580 




6590 




6703 


2S 


1690 




2689 




5861 




6840 




6876 



a In Ref. 
b In Ref. 



SS 

221 



As it can be seen the theoretical results obtained are only estimates for the experimental 
values. The main differences between these theoretical results and the experimental ones, 
or other theoretical calculations, are clearly due to the non-singularity of the potential at 
the origin, and its concavity. That is, the harmonic potential has positive second deriva- 
tive at variance with other proposed potentials that vary more slowly with the distance. 
This causes that when the energy of a state increases, the classical allowed region will be 
smaller for the harmonic potential than for the other ones and states result more localized. 
Other undesirable effects are the constant spacing between the consecutive levels and the 
degeneracy present. The differences between the proposed potential and the former ones 
also cause the so called constituent quarks masses to be here like twice of the usual values 
for the heavy quarks c, b and fifth times for the s quarks. The addition of the non harmonic 
term improved the results for the lower levels and broke the degeneracy present, but for the 
higher excited states the deficiencies remain; and this potential was also unable to fit with 
the flavor independence condition. 
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III. PROPERTIES OF THE RADIAL WAVE-FUNCTION AT THE ORIGIN AND 
RELATED MAGNITUDES 



As it was mentioned in the previous section, the major reason for differences between our 
calculation for the pure harmonic oscillator and previous ones are the non- singularity of the 
potential at the origin, and its concavity. The value of the radial wave function or its first 
non- vanishing derivative at the origin 

*2(o) 



d l R n ,i (r) 



dr l 



(14) 



r=0 



is needed for the evaluation of pseudo-scalar decay constants and production rates through 



heavy- quark fragmentation |31l . 



the ones presented for other potentials in Ref. |33J] for cb mesons. 



32J . In Figure El we compare our results for \R$(0) 
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FIG. 2: Radial wave functions at the origin and related quantities ( R n ({0) (GeV) ) for cb 
mesons. 



Two interesting magnitudes to evaluate, related with the wave functions at the origin, 
are the leptonic widths and the hyperfine splitting. 

Leptonic widths for charmonium and upsilon systems are presented in Figure E3 and 
compared with results in Ref. 24, 



T(V C 



They were obtained by the formula 



l6irN c a 2 e 2 q |# 



Ml 



(15) 



where A^ c = 3 (number of colors), a denotes the fine structure constant, e q denotes the quark 
charge and My is the mass of the vector meson. 



10 r iK2S) 9-57 
^(1S) 9.10 



i.OO 
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1.71 



0.76 



5.81 
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5.26 



2.14 
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V(2S) 
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T(3S) 459 
T(2S) ~ UJ1J - 
T(1S) 



0.55 0-57 



0-52 T(2S) 



Vf 



0.40 0.48 T(3S) 
EQ[33] FU[24] Expt[23] 



FIG. 3: Leptonic widths (KeV). 
Finally, the hyperfine splitting can be obtained through the expression 



M ( 3 5i) - M (%) = Cte 



(16) 



m a m h 

in which the constant is fixed through the hyperfine splitting observed in the charmonium 
family 

M (J/ij) - M (rjc) = lllMeV. 

Results for the cc (J/^ = 3 S l rj c = l S ), be {B* = 3 Si B c = l S ) and bb (T = 3 Si rj = l S ) 
resonances are shown in Figure 0J 
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FIG. 4: Hyperfine splitting for quarkonium ground states (MeV). 
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IV. BARYON SPECTRA 



In the present section we study the baryon spectra using the harmonic potential obtained 
in Section |HI For the three-body system the Hamiltonian has the form 



H 



— V 2 

2m x ri 



1 " V r 2 - 77^V 2 3 + V 12 (n a ) + V 23 (r 23 ) + V 31 (r 31 ) 



(17) 



2rri2 2 2r«3 

In order to separate the C-M motion, we define the Jacobi (~r, R) and C-M {Rcm 
coordinates through 



R 

Rcm 



Hbc 
Hbc 



m b rt + m c rl 



m b + m c 



m a r a +m b r b + m c r' c 
M 



(18) 



where 



Considering equal masses m a 



^bc 
f^a,bc 

M 
m b = 



m b m c 
m b + m c ' 
m a (m b + m c ) 
M 

m a + m b + m c . 



(19) 



m c = m, and Eqs. ()18|) . (fT9*j) the C-M motion is 



separated and the Hamiltonian for the relative motion takes the form 



2/i 



;Wo, 



(20) 



with 



H = 



m a m b m c \ 3 
M 



m 

7!' 



where the rule adopted by Richard |8( (0), was considered. 

At this point we could obtain the baryon spectra directly from Eq. f!20j) noticing that it 
is the sum of two independent harmonic oscillators. But then the energy eigenvalues will 
not be in terms of the natural quantum numbers of the system, and it will not be possible 
a check with experimental or other theoretical results. The same problem is faced with 
the usual 3-dimensional harmonic oscillator, which could be solved in Cartesian coordinates 
as the sum of three independent 1-dimensional harmonic oscillators, but then no relation 
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between the conserved angular momentum and the energy spectrum is obtained and the use 
of spherical coordinates is convenient. 

n 

Then introducing the hyper-spherical coordinates [36| (see Appendix EJ 



r x 


= p cos (x) sin ( 


r ) COS (<f r ) , 


r v 


= p cos (x) sin ( 


r ) sin (lf r ) , 


r z 


= pcos (x) cos 




Rx 


= psin(x)sin( 


Br) COS ((^if ) , 


Ry 


= p sin (%) sin ( 


fl ) sin (ip R ) , 


R z 


= p sin (x) cos ( 





(21) 

the kinetic term stay in a diagonal form and the potential become only dependent of the 
hyper-radio, then the Hamiltonian (J2(J|) takes the form 



H = 

2/i 



1 d ( 5 9 \ 1 

^ pV +3 



9 



p 5 3p V dp) p 2 \ sin 2 (2x) d\ 



sin 2 ( 2X ) £ ) + 



J*(6 r ,<p r ) , L 2 (6 R ,<p R ) 



cos 2 (x) 



+ 



sin (x) 



+^-k + -Wo, 



(22) 



4 v- / 2 

In which J is the angular momentum of the subsystem be and L is the angular momentum 
of particle a respect to the C-M of the two body subsystem be. 

J d_ f . _d_\ 1 d 2 



J 2 



r/ Sill (0 r ) dif r 

f . n d \ 1 d 2 

sm K - 



(23) 



sm {0 R ) 89 R \""~""dO R J sin 2 (Or) dip R 
The Schrodinger equation in this case has also analytical solution (see Appendix El? with 
eigenvectors 



v I'y.a.„„,././,. (£, X, r , <p r , Or, cp R ) = N N Xj ^ '~ £, x L x + 2 (£ 2 ) cos 1+ i ( X ) sin'+s (x) x 



cos (2 X )) YP (Or, <p R ) Yp (0 r , Ifr) , (24) 



where 

A = 2n + j + I. 
The eigenvalues are given by the expression 



(25) 



E 



«(2tf + A + 3) + 2*r 0> 



(26) 
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where N is the number of nodes of the hyper-radial function, A is the grand-angular quantum 
number, and P^-J-l' +5 (cos (2%)) are the Jacobi polynomials. In Figures El and El we compare 

2 

our results with other calculations presented by Richard in Ref 
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FIG. 6: Baryon sss (left) and bbb (right) energies (MeV). 
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To compare with experiments there were only two equal constituent quark baryons at 
hand j^, for the sss system, for which 



m(n-)^ = 1672 MeV, 



m(Dr) Ea!pt = 2250MeV, 
m [Qt) v fi = 1641 MeV, m (^) v fi = 2225 MeV, 



m(f! ) — mifl ) B 

V / Expt \ >Vf 



m ( n )Bxpt 



0.02, 



mifl ) — m[Cl ) 

\ J Expt V l\ 
m ( n ")i3xpt 



0.01. 
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V. CONCLUSIONS 



In the present paper we have studied, within a non-relativistic spin-independent model 
with harmonic confining potential, the spectra and other properties of hadron systems. 

It was found that for mesons, without imposing the flavor independence condition, two 
possible potentials with harmonic confining terms had analytical solutions that give good 
estimates of the experimental values reported for the meson spectra. The better fit was 
obtained for the potential with a term proportional to \ because it has a singularity for 
r — > that improves its behavior in this region. However this potential was found to be 
incompatible with the flavor independence condition and was not considered in the analysis 
that followed. 

For the pure harmonic oscillator the parameters introduced were fixed from the low ly- 
ing levels of heavy quarks systems and imposing the flavor independence condition. The 
calculation of the meson and baryon spectra, and the hyperfine splitting with this potential 
give good estimates of the experimental and other theoretical results; in the case of leptonic 
widths we could say that the results are not good. Although this potential is far from being a 
good approximation for the real inter-quark potential, and the results are not as good as the 
theoretically obtained by other phenomenological models, it has the great advantage that 
allows to obtain analytical solutions for both meson and baryon spectra. That is, reasonable 
theoretical results are obtained without the need of numerical methods and computational 
calculations. The major differences of this potential and the others mentioned in the intro- 
duction are due to its non- singularity at the origin and its concavity, that cause the bad 
results for the obtained leptonic widths and also (for a better fit with experiments) the so 
called constituent masses to be bigger than the usual ones. The baryon spectra was studied 
with the use of the rule (J2J) for the qq potential and for obtaining the analytical solution it 
was necessary to restrict the study to equal constituent quark systems. 

We finalize with the conclusions of A. Martin in Ref. Q| "... if you want to know the 
mass of a particle and if you have a little time (in years!) and little money you should forget 
all your prejudices and use potential models". 
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APPENDIX A: THE HYPER-SPHERICAL COORDINATES AND THE SOLU- 
TION OF THE THREE-BODY PROBLEM WITH A HARMONIC POTENTIAL 



The hyper-spherical coordinates are very useful for dealing with the three-body problem; 
in what follows we make a small review of them. 
The kinetic energy of the Hamiltonian ([20)1 . 

£ = -^(V T 2 + V^, (Al) 

can be written as a Laplacian in a 6-dimensional space, due to the symmetry in the two Jacobi 
vectors. When a change of coordinates is made, the A r -dimensional Laplacian transforms as 

m 



a-ifcE^viJ where H£©> (A2) 

are the metric coefficients. The change to hyper-spherical coordinates is based on the defi- 
nition of the hyper-radius by 



P 



N 



K (A3) 
\ i=i 



and N — 1 angles in a way that (jA3|) is satisfied and the old variables are expressed in terms 
of the new ones by N functions with the form 

Xi = pFi (ftjv-i) , (A4) 

and the Laplacian operator becomes 

p N ^ dp V dp) p 1 

In the expression (jA5|) VIn-i denotes all angles. The explicit form of the angular term 

of the Laplacian operator and its eigenfunctions will depend of the set of angles selected as 

new coordinates and the eigenvalues will be equal to —A (A + N — 2). In this 6-dimensional 

case where the selected coordinates are (|21j) the angular term obtained is 

a 1 9 ( 2 . . d \ 1 1 . . 

sm 2 (2 X )dx\ dxj cos 2 (x) ' v sm 2 ( x ) 

Its eigenfunctions are expressed as a product of orthogonal polynomials in separated 

variables The eigenfunctions corresponding to Jacobi's angles 8 r , tp r , Or and <pn are the 
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well known spherical harmonics. Then the equation obtained for the function of the angle 
X is 







sin 2 (2x) dx 



sin 2 2Xtt- 



1(1 + 1) + 



*V (X) = -A (A + 4) F? ( X ) , (A7) 



dx) sin 2 (x) cos 2 (x). 
for which the solution are the Jacobi polynomials 

(x) = N x , u sin' + 5 ( x ) cos^ + § ( X ) 



1,3 



2 « ' 2 
A-j-i 



(cos (2 X )) 



(A? 



with A = 2n + / + j and N\ij a normalization factor. 

For the harmonic interaction between equal mass particles, the potential is only dependent 
of the hyper-radius, then we can separate variables and the radial equation has the form 



1 d ( 5 d \ A (A + 4) s/3k 



2pp 5 dp \ P dp) 2/ip 2 



Introducing new variables 



and the new function 



we obtain for (|A9J) the equation 



P 



Rn,x (p) = 0. (A9) 



2p (E - |Wo) 



(kV3) 



8p (E - |Wo) 



2 • 



T, 



R 



N,X 



AT, A 



.5 ' 
■ 2 



(A10) 



(AH) 



d 2 



+ 



A(A + 4) + f . 



d£ 2 P. 



N,X 



T N ,x (0 = 0, 



with solutions 



where 



4s + 2 (A + 2) + 2, 



(A12) 



(A13) 



(A14) 



AT, A 



and L^ 2 (£ 2 ) are the generalized Laguerre polynomials. 

Then the complete eigenvectors and the eigenvalues of the problem are given by Eqs. 
and (JUJ). 
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